Relaxation of non-order parameter field in directed Ising systems 
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We investigate the effect of initial conditions on the dynamic exponents of the interacting 
monomer- monomer model with infinitely many absorbing states in one dimension. This model 
exhibits a directed Ising (DI) type transition from an active phase into an absorbing phase. In case 
of the directed percolation universality class, it has been reported that the non-order parameter as 
well as the order parameter exhibits critical fluctuations, relaxing algebraically to its natural value 
with the same scaling exponents. We numerically confirm that this is also valid for the DI univer- 
sality class. We also observe continuously varying dynamic exponents with a linear dependence on 
the non-order parameter initial density. 
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Various kinds of nonequilibrium lattice models exhibit- 
ing absorbing phase transitions have been studied exten- 
sively during last few decades |ffl). Two distinct types 
of absorbing phase transitions have been identified in 
one dimension: the directed percolation (DP) and di- 
rected Ising (DI) universality class. Most models have 
been found to belong to the DP class, which involves 
typically a single absorbing state or multiple absorbing 
states without any symmetry. Dl-type critical behavior 
appears in models with two equivalent absorbing states 
or two equivalent group of absorbing states 



Recently, the question has been addressed whether one 
can construct initial states that affect the entire tempo- 
ral evolution of nonequilibrium systems including asymp- 
totic dynamics. Typical examples are the systems that 
display an absorbing phase transition from an active 
phase into an absorbing phase with infinitely many ab- 
sorbing (IMA) states. As in the ordinary order-disorder 
phase transition models, the order parameter takes a 
nonzero value in the active phase and vanishes in the ab- 
sorbing phase. The IMA states can be characterized by a 
non-order parameter. In the steady state, the non-order 
parameter approaches algebraically to a natural value at 
criticality. 

Nonuniversal dynamic properties have been reported 
in various kinds of DP-type nonequilibrium lattice mod- 
els with IMA states The order-parameter dy- 
namic exponents vary continuously with the initial con- 
ditions characterized by the non-order parameter density. 
Moreover, at criticality, the non-order parameter density 
also exhibits critical fluctuations, relaxing algebraically 
to its natural value with the same dynamic exponents 
0. The critical relaxation of the non-order parameter 
field in DP-type systems with IMA states was recently 
examined from a field theoretical (Langevin equation) 
approach |^,[^] . The evolution equation for the order pa- 
rameter is non-Markovian and includes a temporal mem- 
ory term due to the non-order parameter field. This 
non-Markovian term is known to be responsible for the 



nonuniversal dynamic exponents. 

In this paper, we introduce an interacting monomer- 
monomer model with IMA states (IMA- IMA model) that 
belongs to the DI universality class and investigate the 
effects of the non-order parameter field in one dimension. 
Using numerical simulations, we show that the non-order 
parameter plays the same role in DI systems as well as 
DP systems. 

The IMA-IMM model is an interacting monomer- 
monomer model with two different species of monomers, 
namely A and B. Monomers A and B are selected with 
probability p and 1 — p respectively. Monomers can ad- 
sorb at a randomly selected vacant site with unit proba- 
bility for A and 1 — re for B respectively. An adsorption 
attempt of a monomer is rejected when both sites adja- 
cent to a selected vacant site are occupied or at least one 
adjacent site is occupied with a monomer of the same 
species. A nearest neighbor AB pair reacts and desorbs 
immediately from the lattice. 



A + V 

B + V 
A s + B s 



(l-p)(l-r b ) 



A s 
Bs 



(1) 



Here the subscript s denotes adsorbed particles, V de- 
notes a vacant site where the adsorption attempt is al- 
lowed. 

This model has infinitely many absorbing states. Any 
configuration without a nearest neighbor pair of vacant 
sites is one of the absorbing states. The number of ab- 
sorbing states diverges exponentially with system size L 
as 2 L / 2 . The absorbing states can be divided into two 
equivalent groups; odd-site occupied and even-site occu- 
pied groups. Clearly, these two groups of absorbing states 
have one-to-one correspondence and the Ising (Z%) sym- 
metry in between. So we expect that this model exhibits 
the Dl-type absorbing phase transition. 

We locate the critical line in the r^—p phase diagram by 
dynamic Monte Carlo simulations |n]]. For small values 
of f;,, the system is always absorbing. As rb increases, 
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a window of the active phase appears and divides the 
absorbing phase by A-dominated and B-dominated ab- 
sorbing phases. For example, as p increases along the 
rt = 0.9 line, the system undergoes two continuous phase 
transitions from the i?-dominated absorbing phase into 
the active phase at p c \ = 0.100(3) and finally into the 
A-dominated absorbing phase at p C 2 — 0.5058(7). We 
numerically confirm that both transitions are of the DI 
type as expected jllj . 

The order parameter p is the number density of near- 
est neighbor pairs of vacant sites, while we define the A 
particle density as the non-order parameter pa- First, 
we perform static Monte Carlo simulations to measure 
the natural density of the non-order parameter at both 
criticalities. We measure the A particle density pA{t, L), 
averaged over 2 x 10 3 ~ 5 x 10 4 survived samples for 
system size L = 2 6 ~ 2 10 . pA(t, L) relaxes to the natural 
density p A at = p ,1(00,00) in the thermodynamic limit. 
At p = pd and p — p c2 , we estimate p n A at — 0.101(1) 
and 0.458(2), respectively. The distribution function for 
the A particle density appears to be Gaussian, so the 
mean density is identical to the most probable density. 
We measure the most probable density in the long time 
limit, which turns out to be consistent with the above 
mean natural density within statistical errors. 

Dynamic properties for the non-order parameter 
PA(t, L) can be extracted by studying its temporal de- 
viation from its steady-state value (natural density) as 



Ap A (t,L) = \p A (t,L) 



(2) 



As in the DP case, we assume that ApA(t, L) follows the 
same scaling behavior as the order parameter. Using the 
finite-scaling theory JlJ] in the steady state, Ap A (t, L) 
scales at criticality as 

Ap A {oo,L)~L-V u ±. (3) 

One can also expect the critical short time behavior as 

Ap A (t,oo)~t-^, (4) 

and the characteristic time ta{L) scales as 

T A {L)~mi"±. ( 5 ) 

From static simulations, we estimate the scaling ex- 
ponents at both criticalities. At p = p c i,Pc2, we esti- 
mate 0/v± = 0.46(2), 0.46(3) 0/^ = 0.26(2), 0.25(2), 
and v\\/ux_ = 1.75(5), 1.90(10) (see Fig. 1). As expected, 
these estimations involve rather large statistical and sys- 
tematic errors, especially due to inaccuracy of the natural 
density values. However, these values agree reasonably 
well with the DI values p3[ , which confirm our assump- 
tion that the non-order parameter exhibits the same type 
of critical fluctuations as the order parameter. 



Dynamic exponents for the order parameter in DP sys- 
tems with IMA states are known to depend on initial con- 
ditions characterized by the non-order parameter. To in- 
vestigate this nonuniversal dynamic properties in DI sys- 
tems, we perform dynamic Monte Carlo simulations with 
various initial conditions. We start with a pair of nearest 
neighbor vacant sites in the absorbing background which 
is controlled by the A particle density p A . 

We measure the survival probability P(i)(the probabil- 
ity that the system is still active at time t) and the mean 
number of pair of vacant sites (order parameter) N(t) 
averaged over all samples. At criticality, these quantities 
scale algebraically in long time limit as 14 



P(t) ~ r 



N(t) ~ t". 



(6) 



The dynamic exponents are in general functions of initial 
non-order parameter density; S — S(p A ) and r\ = rj(p A ). 
With initial configurations of the natural non-order pa- 
rameter density, the exponents take the ordinary DI val- 
ues; 5{p n A at ) =S DJ ^ 0.285 and ri{p n A at ) =T] DI ~ 0.00. 

At p — pd with rb — 0.9, we estimate the values of 
5 and rj as p A varies from to 0.35 (see Table I). In 
Fig. 2, we plot the exponent shifts from the DI values 
(5 — Sdi and t] — rjni) versus p\ — p A at '. It shows a lin- 
ear dependence of the exponent shifts on the deviation 
of the non-order parameter from the natural density in 
initial configurations. This linear dependence has been 
also seen in DP systems. 

In summary, we investigated the IMA-IMM model in 
one dimension, which show the Dl-type continuous phase 
transition from an active phase into an absorbing phase 
consisting of infinitely many absorbing states. We show 
that the non-order parameter exhibits critical fluctua- 
tions identical to the order parameter. Dynamic expo- 
nents 5 and r\ depends linearly on the initial non-order 
parameter density and coincide with the ordinary DI val- 
ues only at the natural density. 

This work was supported by an Inha University re- 
search grant. 
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FIG. 1. Temporal dependence of ApA and size depen- 
dence of ApA and ta at p — p c \. The straight lines are of 
slope 0.26(=/3/i/||), 0.46(=/?/V±), and 1.80(=V||/i/±). 
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FIG. 2. Initial configuration dependence of the exponents 
S, r\ at p = pd. It shows a linear dependence on the initial 
non-order parameter density. 



TABLE I. Initial configuration dependence of the dynamic 
exponents at p c — 0.1. 



PA 


8 


V 


0.00 


0.352(5) 


-0.065(10) 


0.05 


0.322(5) 


-0.033(7) 


0.10 


0.293(4) 


-0.006(8) 


0.15 


0.258(8) 


0.031(9) 


0.20 


0.223(4) 


0.060(7) 


0.25 


0.193(5) 


0.090(7) 


0.30 


0.155(3) 


0.124(6) 


0.35 


0.122(4) 


0.155(8) 
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